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Recent efforts in cosmic ray (CR) confinement and transport theory are discussed. Three 
problems are addressed as being crucial for understanding the present day observations and 
their possible telltale signs of the CR origin. The first problem concerns CR behavior right 
after then - release from a source, such as a supernova remnant (SNR). At this phase the 
CRs are confined near the source by self-emitted Alfven waves. The second is the problem 
of diffusive propagation of CRs through the turbulent ISM. This is a seemingly straight¬ 
forward and long-resolved problem, but it remains controversial and reveals paradoxes. A 
resolution based on the Chapman-Enskog asymptotic CR transport analysis, that also in¬ 
cludes magnetic focusing, is suggested. The third problem is about a puzzling sharp (~ 10°) 
anisotropies in the CR arrival directions that might bear on important clues of their transport 
between the source and observer. The overarching goal is to improve our understanding 
of all aspects of the CR’s source escape and ensuing propagation through the galaxy to the 
level at which their sources can be identified observationally. 


I. INTRODUCTION 

Cosmic rays (CR) have been discovered more than a century ago but the problem of their origin 
is still with us today. The fundamental obstacle to identification of their possible sources, such 
as the supernova remnant (SNR) shocks, is the CR ’black-box’ propagation through the chaotic 
magnetic field of the galaxy. With a possible exception for the highest energy CR (> 10 19 eV, 
whose origin is almost certainly extragalactic [ 1 ,0) most of the CRs arrive from random directions 
saying nothing about the locale of their sources. The more surprising is a sharp (~ 10°) CR 
anisotropy discovered by Milagro [|3] with interesting ramifications due to IceCube, ARGO-YBJ 
0], HAWC [0 and some other instruments. Had it been created in the source, it would have 
been completely erased en route to the Earth. Therefore, this CR feature is likely to be an imprint 
of their interaction with the ISM (interstellar medium) or its local environment (heliosphere and 
surroundings). We will discuss this later. More natural is to start with the CR transport right at 
their birth place. 

SNRs are widely regarded as the most probable source of the bulk of the CRs [2]. The joint 
analysis of the broad band observations of the SNRs and CR background spectra at the Earth 
should provide the ultimate evidence for this hypothesis. The analysis faces multiple problems. 
First, the accelerated CRs manifest themselves in SNRs only in form of secondary emission, which 
is usually difficult to interpret. For example, the super TeV-photons, carefully counted by the at¬ 
mospheric Cerenkov telescopes, to testify for the accelerated protons colliding with the ambient 
gas L6, 0 , can easily be confused with the inverse Compton (IC) photons up-scattered by accel¬ 
erated electrons. If this is the case, not much weight can be added to the argument for the CRs 
origin in the SNRs, as electrons comprise only a small fraction (~ 1 — 2%) of the CR spectrum. 
Similarly, in the GeV energy band the emission may come from electron Bremsstrahlung. The key 
in both cases, however, is a dense gas in the SNR surroundings, often present in form of adjacent 
molecular clouds (MC). They provide a target for the pp reactions with accelerated protons. Pho- 
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tons, produced in this reaction should thus come from MCs, illuminated by the CR protons that 
have, in turn, escaped from the source. By contrast, the IC electron emission should come from 
the entire volume they fill, as the low-energy photons (such as CMB) are present everywhere. To 
use this simple but powerful diagnostic tool for identification of the source of emission detected 
by modern ground based instruments and space observatories (such as the Cerenkov teles cop es 


HESS, VERITAS, MAGIC and Fermi-LAT, PAMELA, Agile spacecraft observatories |[8l4l4j]). 
we must have an accurate understanding of the CR propagation from their sources to the adjacent 
MCs. 

The next problem is the subsequent interaction of the CRs with the MC, as their confinement 
inside the cloud is generally deteriorated due to the collisional damping of the Alfven waves, that 
otherwise would prevent CRs from spreading further rapidly. In addition, this interaction reveals 
important clues as to how the spectrum of CRs, illuminating the MC is different from that in the 
source, most importantly, in form of spectral breaks. This aspect of the CR interaction with MC 
and their visibility in the gamma-ray band has been discussed in one of the earlier APS Plasma 
Physics meetings [15fl. Here we will focus on the ensuing propagation of the CR to the Earth and 
their spectral features that they can acquire both in the source and on the way to us. 

Measurements of the CR background spectrum have also advanced significantly. Much 
progress has been made in isolating different elements in it. One of the most striking results 
was the ps 0.1 difference between the rigidity (momentum to charge ratio) spectral indices of pro¬ 
tons and He 2+ ions. Such deviations have been apparent for some time, e.g., [ f6] but the Pamela 
spacecraft observatory measured it with a three-digit accuracy in the 100 GeV energy band [12], 
which posed a strong challenge to the CR acceleration and propagation models. Indeed, the ultra- 
relativistic parts of the rigidity spectra must be identical, if protons and He 2+ ions are accelerated 
and transported via electromagnetic interactions under identical conditions. Since He 2+ has a 
0.1-flatter spectrum, the difference may be due to its spallation, biased for lower energies [17]. 
However, such scenario would probably require stretching the model parameters too much [18]. 
Other interesting possibilities discussed in the literature include the contribution from multiple 
SNRs of different types with somewhat different CR spectra ltl2i.ll9l-211 and variable /;/Hc 2 mix 
along the shock path [122,[23b. A ’plasma physics’ solution that targets the nonrelativistic phase of 
acceleration of both species, where the argument of equal rigidity spectra is irrelevant, was sug¬ 
gested in Rcf.[24j], This explanation is advantageous according to Occam’s razor, as it relies on the 
collisionless shock intrinsic properties and does not require any of the above special conditions. 

Apart from the elemental composition, other spectral signatures, such as the spectral hardening 
above E ~ 200 GeV, have been studied EO, 12511 and provided important clues for the energy 
dependent CR transport. Stochasticity of CR sources and inhomogeneity of transport through the 
galaxy are now also included in the models [26j]. These are important for understanding the large 
scale CR anisotropy. The most puzzling aspect of the anisotropy in the CR arrival directions is, 
in my view, the sharp anisotropy or the so-called Milagro ’hot spots’ which I address later in this 
brief review. 

The remainder of the paper is organized as follows. In SecjlQthe confinement of CRs released 
from the source is addressed. In SecOH] an equation describing diffusive propagation including 
a hyperdiffusive term is presented and its relation to the so-called “telegrapher” term in the CR 
transport equation is clarified. In Sec[TV]possible mechanisms for building a sharp CR anisotropy 
during their propagation from the source are discussed. 
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II. SELF-CONFINEMENT OF CRS AROUND SNRS 


It is widely believed that CRs are accelerated in SNR shocks by the diffusive mechanism 
(DSA). The backbone of the DSA is a self-confinement of accelerated particles supported by their 
scattering off magnetic irregularities that particles drive by themselves while streaming ahead of 
the shock. Logically, this process should also control the ensuing propagation (escape) of CRs, 
at least until their density drops below the wave instability threshold. At the same time, no con¬ 
sensus has been reached so far as to how CRs escape the accelerator. The dividing lines seem to 
run across the following issues: (i) does the escape occur isotropically or along the local magnetic 
field? (ii) does the scattering by the background MHD turbulence control the CR propagation 
alone or self-excited waves need to be included? (iii) are CRs, that escape SNR, peaked at the 
highest energy or lower energy CRs escape as well? [27 429(1 [301 B ill 02211 [321. 

Adhering to the self-confinement idea, we consider the model that explicitly includes the self- 
excited waves. Moreover, in the regions where magnetic perturbations are weak, i.e. 8B 2 /B 2 <C 
1, a field aligned CR transport dominates, as the perpendicular diffusion is suppressed, Kj_ ~ 
( 8B/B ) 2 Kg <C Kj| — Kb (8B/B)~ 2 . Here Kb is the Bohm diffusion coefficient Kb — cr g /3 with r g 
being the particle gyroradius. Taking into account the condition (5R 2 /R 2 ) ism -C 1, such regime 
is inevitable outside the source where 8B/B < 1, as well as at later times of CR propagation, 
when they are spread over a large volume and the waves are driven weakly. Moreover, the self¬ 
confinement of CRs propagating away from the accelerator, as described below, is the continuation 
of physically the same process long believed to be at work inside the accelerator, as first suggested 
by Bell [33j]. From a mathematical standpoint our treatment below generalizes Bell’s steady state 
solution, obtained in the shock frame, to the time dependent solution for the CR cloud expanding 
further out. This being said, we use the following equations that describe the CR diffusion and 
wave generation self-consistently [34] 


d , d Kb dP c r 

bM p)= SzI dz 

(1) 

d dPcR 

a/ = " c * * - n ' 

(2) 


where Ca is the Alfven velocity. The dimensionless CR partial pressure / j cr is used instead of 
their distribution function f(p,t ): 


„ 471 2 4 f 

PCR = —~7^ V P f , 


(3) 


3 pC A 

where v and p are the CR speed and momentum, and p- the plasma density. The total CR pressure 
is normalized to d In p, similarly to the wave energy density /: 


(8B 2 ) B% f . . B 2 f 

= l(k)d]nk=-± I 

8 K 8 K J 8 K J 


(p)d\np 


The last relation implies a simplified wave-particle resonance condition, kr g ( p ) = const ~ 1. Most 
of the works on CR self-confinement (see [35] for a review) use equations largely similar to eqs.dH- 
[2]) but different assumptions are made regarding geometry of particle escape from the source, 
the character and strength of the wave damping F, and the role of quasilinear wave saturation. 
A reasonable choice of the damping mechanism is the Goldreich-Shridhar (GS) MHD cascade 
[36], which seems to be appropriate in / < 1 regime [37-39]. The damping rate in this case is 
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r = CAy/kJl , where / is the outer scale of turbulence, which may be as large as 100 pc (see, 
however, SecUVl). As F does not depend on I and can be considered also as coordinate independent, 
it allows the following (’quasilinear’) integral of the system given by eqs.([B and ©: 


Pcr (z,t) =Pcro(z) 


kb d 7(z,f) 
C A dz / 0 (z) 


(4) 


Here P CR0 (z) and 7 0 (z) are the initial distributions of the CR partial pressure and the wave energy 
density (see [34] for more general treatment). Substituting Pcr in eq.© we arrive at the following 
diffusion equation for / 


dl _ d Kg dl dPcm 

dt~ dz I dz A dz 


Outside of the region where Pcro f 0, the last term on the r.h.s. is absent, while the second term 
may be eliminated by replacing I exp (Ft) —> /, exp(Ft)dt —y t. However, if T is taken in a 
GS-form, it is fairly small due to the factor \Jr g jl <C 1. We may simply neglect it. The solution 
for / and Pqr ( z,t ) may be found in an implicit form (see [134)] for details). However, there exists a 
very accurate convenient interpolation formula that can be represented as follows 


Pcr 


2kb(p) 

c) /2 V5 


C 5/3 + (Onl ) 5/6 


3/ Vf 2 /40,SM 


(5) 


where a is the size of the initial CR cloud, £ = z/\/C. A at. Z) NL = F (n) • AsM ex P (— n), with n 
being a normalized integrated pressure 


n = — f p CR dz 
k b J 
0 


The function F, behaves as follows F (n) ~ 2e ~ 5.4, for n 1 and F (n) ~ 27ZTI 2 , for n< 1. 
Here Asm is a normalized background diffusivity Dism = Kb/oCaIism- 

To summarize these results, the self-regulated CR escape from a source is characterized by their 
distribution (partial pressure) comprising the following three zones: (i) a quasi-plateau (core) at 
small z/y/t < i/Av l of the height ~ (DnlO 1,2 • It i s elevated by a factor ~ n 1 exp (n/2) » 1, 
compared to the test particle solution because of the strong quasi-linear suppression of the CR dif¬ 
fusion coefficient with respect to its background (test particle) value Asm: D\jl rs./ Asm exp (-n) 
(ii) next to the core, where a/AJl <z/\/i < A Asm, the profile is scale invariant, Pcr 1 /z. 
The CR distribution in this “pedestal” region is fully self-regulated and independent of n and 
Asm f° r n 1, (iii) the tail of the distribution at z/\ft > v / Asm is similar in shape to the 
test particle solution in ID but it saturates with n> 1, so that the CR partial pressure is 
°c (Asm/) l/2 cxp ( — z 2 /4D|sm/), independent of the strength of the CR source n, in contrast 
to the test-particle regime in which it scales as °c n, (n < 1). Because of the CR diffusivity re¬ 
duction, the half-life of the CR cloud is increased and its width is decreased, compared to the test 
particle solution. Depending on the functions n (p) and Asm (p), the resulting CR spectrum gen¬ 
erally develops a spectral break for the fixed values of z and t at the CR momentum p determined 
by the following relation: z 2 /t ~ Ail (p) ~ Asm cx P (—n). 
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III. DIFFUSIVE AND HYPERDIFFUSIVE CR TRANSPORT WITH MAGNETIC FOCUSING 


Propagating away from their sources, CRs are pitch-angle scattered on weak ISM magnetic 
irregularities. A seemingly straightforward reduction of kinetic CR description to their spatial 
transport leads to a diffusive approximation which has the well-known defect of causality viola¬ 
tion. There have been attempts at an alternative approach based on the “telegrapher” equation. 
However, its derivations often lack rigor and transparency and had not been performed to the re¬ 
quired (as we show below, fourth order) accuracy. The problem can be formulated very plainly: 
How to describe CR transport by only their isotropic component, when the anisotropic one is 
suppressed by the frequent scattering? 

The angular distribution of CRs is described by the function f(fj,,t,z) 100]: 


df 

dt 


d_ 

dfl 




d£ 

dfl 



( 6 ) 


Here £ — v/lv is the small parameter of the problem, with v being the particle velocity, /- char¬ 
acteristic scale and v- scattering frequency. The dimensionless magnetic mirror inverse scale 
a = —B l dB/dz, z points in the local field direction and is measured in the units of /, time - in 
v _1 , /i is the cosine of the pitch angle, while D (j.i) ~ 1 depends on the spectrum of magnetic 
fluctuations. The isotropic reduction scheme requires a multi-time asymptotic (Chapman-Enskog) 
expansion. So, we introduce a set of formally independent time variables t -A to, so that 


d d d 2 d 

3- — 3-b £3-b £ 3— • ■ ■ 

dt dto dt\ dt2 

which leads to the following hierarchy of equations 


(7) 




= ~B 


dfn 


dz 


fo 


A) 


dfn 


d/1 


I 

k=\ 


dfn-k 

du 


( 8 ) 


= -S?b- 1 [/] (t 0 ,...,t n -,iJ.,z) 


where f — fo + £f \ + £ 2 /2 + ■ • ■ and the conditions f n< 0 = 0 are implied. Using the above expan¬ 
sion, one may obtain an equation for the isotropic part /o = (/) = ( 1/2) f fdji to arbitrary order in 
£. By construction, in no order of approximation will higher time derivatives emerge, as was obvi¬ 
ously devised in the Chapman-Enskog method. We terminate this process at the fourth order, £ 4 . 
This is the lowest approximation required to clarify the origin of the telegrapher equation. Higher 
order terms can in principle be calculated at the expense of involved algebra but such calculations 
would be of no avail. So, our result is as follows [42] 


dfo 

dt 


-d'{K-ed''(iiW 2 ) 


(d'lf{W 2 (K-U>))-Ud z 


»(1-M 2 ) / 


dfo 

dz 


(9) 


where d' z — d z + o and d z — d z + a/2. The coefficients are defined as dW / d/l — 1 /D, (W) — 0, 
U'{il) = dU/dil = (l — B 2 ) /D, U (—1) = 0, and K = (1/2) U (1). 
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1. Producing the telegrapher term 


Within the employed Chapman-Enskog expansion, the equation for /o remains evolutionary in 
all orders of £, so no telegrapher term appears. Such term is usually obtained either without clear 
ordering, e.g. [43], using specific D (/i), e.g. [44], or by truncation of eigenfunction expansion, 
where the discarded terms may be of the same order in small parameter as those retained, e.g., 
S0. In most of these treatments care has not been exercised to eliminate the short time scales 
which are irrelevant to the long-time evolution of the isotropic part of the CR distribution, sought 
by these reduction schemes. Instead, they retain the second time derivative which changes the 
type of the resulting transport equation to hyperbolic. As we show below, the second order time 
derivative term can be recovered from the Chapman-Enskog expansion. 

To resolve the above controversy, we simplify eq.© by removing terms unimportant for the 
controversy. First, we may set a = 0 and assume the scattering symmetry, D{—fi) — D(fi), to 
remove the term ~ £ 3 , as the d 3 /dz 3 term is not included in the telegrapher equation derived for 
magnetic focusing by, e.g., [43]. Using these simplifications and the slow time T = £ 2 //4, eq.© 
rewrites: 


dfo „ d 2 / 0 

dz 2 


dT 


= K 


-to?* 


dz 4 


where K is the hyper-diffusion coefficient 





2 D(\-p 2 ) 

To the same order in £ <C 1, the last equation can be rewritten as follows 

dfo d 2 f 0 _d 2 fo 


( 10 ) 


(ID 


= K- 


— T- 


dT “ dz 2 ' dT 2 ’ (12) 

where T = e 2 K/k 2 . This equation has indeed the form of a telegrapher equation. However, the 
comparison of eq.dTOl) with, e.g., the telegrapher eq.(15) in Ref. 04311 shows that, the coefficient t 
in eq. (fl2l) is substantially different. The reason is that the equation of Ref.[43] has been obtained 
by a formal iteration not accounting for all the fourth order terms, the telegrapher term actually 
originates from. Note that [44] give an expression for r which is probably consistent with the 
result above, but their t was obtained for a power-law Z)(/i), so the comparison would take extra 
steps, not worth doing in this short review. More importantly, the telegrapher term in eq. (fT2l) has a 
small parameter (at highest derivative). The role of such terms is known from the boundary layer 
problems. They become crucial near and inside the boundary layer, thus determining its structure 
and scale. In the context of the telegrapher equation, the boundary layer translates into the initial 
relaxation phase of the CR distribution. This relaxation is associated with the small scale CR 
anisotropy in f n which quickly decays. It should be noted that if a simplified collision term (BGK, 
or t- approximation) is used instead of the pitch-angle diffusion in eq.©, the telegrapher equation 
can be accurately derived with no recourse to hyperdiffusion [47]. 

To conclude this section, by comparison with the telegrapher equation, the classic Chapman- 
Enskog is a considerably more suitable and flexible tool to describe the long-time CR propaga¬ 
tion, although the telegrapher version (with corrected transport coefficient) may still be useful for 
studying the magnetically focused CR transport, e.g. [48]. Efforts on improving the CR diffusion 
models, where their drawbacks are important, need to address the lower level transport, including 
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anisotropic component of the CR distribution, directly. Recent work can be found in, e.g., [49Q 
and in the next section. Splitting_the particle distribution in scattered and unscattered categories is 
another useful approach, e.g., [501, 51]. However, when the diffusive treatment is well within the 
method’s validity range (weakly anisotropic spatially smooth CR distributions) neither the teleg¬ 
rapher term nor hyperdiffusivity are essential to the CR transport (see 0 for more discussion). 


IV. SMALL-SCALE CR ANISOTROPY 


CR acceleration (e.g. DSA) and propagation models, as discussed in the two preceding sec¬ 
tions, typically predict only a large scale, dipolar anisotropy. It would emerge as a small /i ^ /i 
correction to /o>/i, produced by localized sources, and can be easily obtained within the treat¬ 
ment outlined in SecjllT] The same is true for the CR self-confinement problem considered in 
SecJH if the small anisotropic correction is taken into account. Now that we expect the CR propa¬ 
gation in the essentially stochastic magnetic fields to be largely ergodic, there is no obvious reason 
for a significantly sharper than the dipolar anisotropy. Yet observations show that narrow (~ 10°) 
CR beams do exist [3-[5]. They shed new lights on the CR propagation from, and even their 
acceleration in, putative sources and need to be understood. 

A number of scenarios have been suggested to explain the tightly collimated beams. They in¬ 
clude magnetic nozzle focusing [52i| . propagation effects from local SNR [53|], acceleration in the 
heliotail [54, 55] and heliosheath propagation effects [56]. Although being plausible in principle, 
those explanations impose significant constraints on the relevant parameters and processes. For 
example, the magnetic mirror ratio must be rather strong to produce ~ 10° anisotropy, and quite a 
strong magnetic field in the heliotail is required to confine 10 TeV protons and make the proposed 
acceleration mechanism work efficiently. Conceptually different scenarios 00 essentially 
attempt at generating small-scale anisotropy out of the large-scale one by exploiting aspects of 
interactions between the CRs and MHD turbulence in the ISM. At the first glance, precisely the 
opposite should occur and the task is clearly of a kind of “squeezing blood out of stone”. From 
a purely mathematical perspective, using certain properties of the particle propagator, these mod¬ 
els produce multipoles out of the dipolar component over a long distance (up to a few lOOpc) of 
particle propagation. At this point, however, the approaches deviate strongly from one another. 

In Ref. [590, an interesting technique is employed to generate higher multipoles from the avail¬ 
able dipole by using the Liouville’s theorem. It is not clear, however, whether the introduction of a 
simplified collision term in a BGK-form, is justified for the treatment of the small-scale anisotropy. 
The preferred collision operator is the differential one which is much more efficient at smoothing 
small-scale anisotropies (see, e.g., previous section). Attacking the same problem from a different 
angle, the authors of Ref. [58]] rightly state that, although the scattering fields are random, we do not 
really need to perform an ensemble average, as the current MHD turbulence is static for the limited 
time observations and fast CRs. There are at least two tests to propose for this explanation. First, 
as this is actually a magnetic lensing effect with a very long particle path (L r g ), small variations 
of magnetic configuration may produce significant changes in arrival directions of narrow beams. 
Indeed the relevant scale of the turbulent field is r g , so the time scale is T ~ r g j(V a + Uhs ) with 
the Alfven and the heliosphere velocities in denominator. The median Milagro energy is ~ ITeV, 
so for Va + Uhs —50km/s and B — 4/iG one obtains t < lOyrs. This may be close enough to the 
time difference between Milagro and ARGO/HAWC more recent observations. And yes, changes 
are being observed but they are not quite significant and HAWC is not fully operational yet, so 
more observations are required and they are underway [U]. The second test has, in fact, already 
been performed by the authors of Ref.[58]. Since CRs interact with the static magnetic fields, their 
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dynamics may be regarded as almost ergodic (strong orbit mixing) on every isoenergetic surface 
in phase space. Small deviations from ergodicity are responsible for the hot spots in arrival di¬ 
rections. Moving from one energy surface to the next by A E ~ E should strongly decorrelate the 
spots, since A r g ~ r g for them. This is indeed observed in simulations carried out in Ref.[58]. The 
upcoming improvements in the energy spectra measurements [5] should substantiate such tests 
quantitatively and help to discriminate between different mechanisms. 

The approach of Ref.[57] is also based on the CR interaction with the ISM turbulence, but 
includes ensemble averaging, thus removing the above concerns with the short time variability 
(except for the possible heliospheric variations [60]). The beam direction is assumed to be along 
the local large scale magnetic field (// oc r g ), to minimize the curvature and gradient drifts, that 
would otherwise evacuate particles from the magnetic tube connecting observer with the source, 
since the drifts increase with the pitch angle. The following assumptions are made to obtain the 
beam collimation: (i) large scale anisotropic distribution of CRs (generated, for example, by a 
nearby accelerator, such as a SNR, magnetically connected with the Earth) and (ii) Goldreich- 
Shridhar [61J] (GS) cascade of Alfvenic turbulence originating from a specific scale /, which is the 
longest scale relevant to the wave-particle interactions. 

It is found that the CR distribution develops a characteristic angular shape consisting of a large 
scale anisotropic part (first eigenfunction of the pitch-angle scattering operator) superposed by a 
beam, sharply focused in the momentum space along the local field. The large scale anisotropy 
carries the momentum dependence of the source. The following four quantities are tightly con¬ 
strained by the turbulence scale /: (1) the beam angular width, that increases with momentum 
as oc yfp (2) its fractional excess (with respect to the large scale anisotropic component), that in¬ 
creases as oc (3) the maximum momentum, beyond which the beam is destroyed via instability, 
Pmax- If the large scale anisotropy originates from a nearby source, magnetically connected with 
the Earth, the model predicts (4) the range of possible distances to this source, Is ~ 100 — 200pc. 
If such source is absent, this range corresponds to the beam collimation length, also a few lOOpc, 
with the large scale anisotropy originating from the smooth omnigalactic CR gradient. This scale 
is consistent with the beam collimation length, obtained numerically in [58]. 

If the turbulence outer scale l is considered unknown, it can be inferred from any of the first 
three quantities (1-3) as measured by MILAGRO. All the three quantities consistently imply the 
same scale l ~1 pc. The calculated beam maximum momentum encouragingly agrees with that 
measured by MILAGRO (p max ~ 10 TeV/c). The theoretical value for the angular width of the 
beam is found to be A$ ~ 4 y/e, where e — r g (p) /1 < Cl. The beam fractional excess related to the 
large scale anisotropic part of the CR distribution is ~ 50e. Both quantities also match the Milagro 
results near its median energy, that is for £ ~ 1 — 2 TeV. So, the beam has a momentum scaling 
that is one power shallower than the CR carrier, it is drawn from. One interesting conjecture from 
the l ~1 pc requirement is that the proton ’knee’ at ~ 3 PeV and the beam are of the same origin, 
as these particles may provide the required outer scale for the MHD turbulence, r g ~ /. Another 
possibility is to employ the spiral-arm 1-pc value for /, as suggested in [62]. 
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